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SUMMARY 


In this paper the author discusses the problem of optimization of shallow 
frame structures which involve a coupling of axial and bending responses. A 
shallow arch of a given shape and of given weight is optimized such that its limit 
point load is maximized. The cross-sectional area, A(x) and the moment of inertia, 
I(x) of the arch obey the relationship I(x) = p [A(x)] n , n = 1,2,3 and p is a 
specified constant. Analysis of the arch for its limit point calculation involves 
a geometric nonlinear analysis which is performed using a corotational formulation. 

The optimization is carried out using a second— order projected Lagrangian 
algorithm and the sensitivity derivatives of the critical load parameter with 
respect to the areas of the finite elements of the arch are calculated using 
implicit dif f erentation. Results are presented for an arch of a specified rise to 
span ratio under two different loadings and the limitations of the approach for the 
intermediate rise arches are addressed. 


INTRODUCTION 


With the advent of highly flexible large space structures the nonlinearity of 
response of such structures plays a dominant role in the control of such structures. 
Naturally, optimization of structures in nonlinear response is gaining prominence. 
This paper addresses the issue of optimizing shallow frame structures in nonlinear 
response involving a coupling of axial and bending actions. The objective is to 
optimize a shallow arch of a given shape and given weight such that its limit point 
load is maximized. Besides having to perform a nonlinear analysis in calculating 
the limit point load an issue of even greater concern is that of calculating the 
sensitivity derivatives of the critical load parameter with respect to the design 
variables, namely the cross-sectional areas of the elements of the discretized model 
of the arch. Two approaches are available for the calculation of sensitivity 
derivatives: the direct and the adjoint approach [1]. In general, the adjoint 
approach is preferred for problems involving nonlinear response [2] - [A], The 
popularity of the adjoint approach stems from the fact that the differential 
equations governing the adjoint variable are linear even though the corresponding 
equilibrium equations in terms of the true displacement variables are nonlinear. 
But to date the author is unaware of the use of the adjoint approach for problems 
involving limit point instability. The present work outlines a direct approach 
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similar to that used by the author and his co-worker in the case of shallow space 
trusses [5]. In this approach the sensitivity derivatives of the critical load 
parameter are obtained through an implicit differentiation of the nonlinear 
equilibrium equations as explained below. The present discussion is restricted to 
finite-element models of shallow arches whose cross sections obey the relationship 

I(x) = p [A(x)] n , n = 1,2,3, P = specified constant (1) 


SENSITIVITY DERIVATIVES OF THE CRITICAL LOAD PARAMETER 


Consider a shallow arch under a given distribution of loading. Assume that 
i cr is the smallest value of the parameter by which the given distribution of 
loading must be scaled in order to produce instability of the arch. The parameter 
i is then defined by the solution of the following system of equations of a finite- 
element model of the arch. 
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where n denotes the total potential energy of the model undergoing finite 
displacements and qj , 1 = 1,2 ...N denote the generalized nodal displacements of 

the model. The load parameter X occurs implicitly in Eqs. (2) and (3). Assume 
that A v for k = 1,2... m are the m design variables, which for the arch are the 
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cross-sectional areas of the finite elements. To obtain 7^ we proceed as 
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follows : 


Rewrite Eqs. (2) and (3) as 
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i , H = 1 , 2 ... N 
k = 1 , 2 ... m 


An implicit differentiation of Eqs. (4) and (5) with respect to A k leads to 
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for every k * 1,2 ...m. These equations may be written symbolically as 



where H is the Hessian matrix of the potential energy of the finite-element model 
of the arch, F is the given vector of nodal forces, and G is the row matrix of 
derivatives of the determinant g of the Hessian matrix with respect to nodal 

displacements. Equations (7) assume that is equal to zero since for constant 
directional loading parameter A does not occur explicity in the stability 
criterion. The elements of G can be evaluated by using the formula 
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where 
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is the matrix obtained by differentiating each element of the deter- 


minant of H with respect to a typical component q } . With this, the sensitivity 
3 A 

derivative can be calculated by the solution of Eqs. (7) at a given (q,A). 
k 


Incidently , Eqs. (7) apply everywhere along the loading path including the limit 
point. It is only at a bifurcation point that the determinant of H is not differ- 
entiable. For very shallow arches instability typically occurs through snap- 
through and hence Eqs. (7) clearly apply. 
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SPECIALIZATION TO A FRAME ELEMENT MODEL OF THE ARCH 


We illustrate the derivation of the sensitivity derivatives for a finite 
element model with 3-D frame elements. For its kinematic description the frame 
element uses the co-rotational formulation as outlined in great detail in reference 
[6]. According to this formulation, which permits large rigid body motion of the 
element, the total motion is decomposed into a rigid body component and a strain- 
producing component. For an element p-q of length L, the displacements of the end 
q relative to the end p in the body fixed axes can be shown to be 
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Uj , Vj , and Wj (i = p or q) denote the global displacements of the nodes and 
matrix [T] p is [6]. 
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c. « cos a., s. = sin a. for i = x,y and z. Angles <|> ,<J> and <b are the initial 
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orientation angles of the frame element and the angles 0 , 0 and 0 are the 
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rigid-body rotations of the end p. In deriving Eqs. (10) Euler angle 
transformation is implied with the order of the rotations being a^, a y and a^. 

Similarly, with the restriction of small relative rotation within the element, 
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Assuming the relative axial and transverse displacements to be linear and 
cubic, respectively , the strain energy of the (p-q)th or the e-th element, e = 
1,2...m can be shown to be [6]. 
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All the expressions 


for the evaluation of matrices in Eqs. (7) are now available and, in principle, can 
be evaluated even though the algebra may be rather tedious. The above expressions, 
especially the [T] p matrix, can be simplified using the assumption of small rigid 
body motions within a load step. 

Indeed, Updated Lagrangian formulation for the kinematic description may have 
simplified matters quite a bit especially if the expressions are linearized within 
a load step but the above expressions using the co-rotational formulation permit 
truly large displacements and with an highly efficient algorithm for the solution 
of nonlinear equations like for instance the BFGS algorithm [7], it can permit 
relatively large load steps resulting in a fewer number of load steps to attain a 
given load level. 


CONSTRAINED OPTIMIZATION 


The optimization problem consists of maximizing the critical value of the load 
parameter A subject to the constraint on total volume of the structure and side 
constraints on member sizes. Although it is perfectly permissible to pose the 
problem as 


min ( f (A) = - A) (15) 
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experience suggests the following well-posed problem 
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where A cr is located by incrementing the load parameter and locating its level at 
which the determinant of the Hessian vanishes. This can be done by monitoring 
either the determinant or the inertia of the eigenvalues of the Hessian matrix H. 
Once an interval is located where the critical point is supposed to lie its 
exact location is determined by a root-finding technique. With Eqs. (15)— (19) 
there is no guarantee that the lowest value of A that satisfies Eqs. (17) will 
always be found. 

The problem as posed by Eqs. (20)-(22) is solved by using Powell's variable 
metric algorithm for constrained optimization (VMCON) [8]. The required gradient 
of the Lagrangian function corresponding to Eqs. (20)-(22) involves the gradient 
of the load parameter which is calculated using the expressions derived in the 
previous section. 


DISCUSSION OF RESULTS AND CONCLUSIONS 


The first step was to validate the accuracy of the sensitivity derivatives. 
This validation was performed by comparing the analytically calculated derivatives 
using the expressions ( 6 ) — ( 8 ) with those calculated using central differences. 
Since no previous studies exist that address the problem being discussed herein, it 
was essential to generate a basis for commparison. Such a basis was provided by 
designs that corresponded to maximum potential energy of the nonlinear deformations. 
Even though previous studies on shallow trusses [5] have confirmed the 
non-optimality of such designs they are relatively easy to generate and provide a 
basis for comparison with truly optimum designs. 

It can be easily verified that designs which correspond to maximum potential 
energy satisfy the condition 

(U s + n U b ) 
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g to 

where U , U and V are the strain energy due to stretching, the strain energy due 
to bending, and the volume of the eth element, respectively. 


Relations (23) can be easily met by a recurrence procedure that evolves 
design for the (r+1)st iteration from that of the rth iteration according to 
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a is a constant such that 
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and p is a suitable exponent usually chosen to be equal to 1/2. Several designs 
for a concentrated load at the crown and a uniformly distributed vertical loading 
were generated for n = 1,2,3 using the mathematical programming procedure, VMCON 
and the recurrence relations (23)- (24). Table 1 provides a comparison of these 
designs. Differences between the two designs are indeed drastic especially for n = 
3. A curious phenomenon was observed during the recurrence procedure namely, that 
several non-converged intermediate designs had higher critical (limit) loads than 
the final fully coverged designs with a uniform specific energy density 
distribution. This is to be expected since the fully converged designs are 
non-optimal. Table 2 provides the material distributions in terms of the 
non-dimensional areas of the five frame elements used to model half the arch. 

tt x 

An attempt to optimize a five element arch model with y(x) = 5 sin failed 

for n = 2,3 because no limit point load could be determined. This is not 
surprising since for very low rise to span ratios the arch is likely to behave more 
like a flexible nonlinear beam with no susceptibility to snap-through. Likewise 
for arches with high rise to span ratios instability occurs by bifurcation at load 
levels far below their limit points and hence the problem belongs to the class of 
linear eigenvalue problems. For arches with intermediate rise to span ratios the 
type of instability can change from the initial limit point to a bifurcation type 
at convergence. In fact the two points may coincide during optimization at which 
point the critical load parameter is no longer differentiable with respect to the 
design variables. Recourse must be then made to techniques of nondif f erentiable 
optimization [8]. 
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Table 1 . Comparison of Designs for Different Loadings on an Arch 


y(x) = a sin — ; a = 10, L = 100 
L 


Type 

of Design 

Ucr^OPT / (*cr) unif* 

Concentrated Load at the 
Crown 

Uniformly Distributed Vertical 
Loading 

n = 1 

n = 2 

CO 

II 

C 

a 

ii 

n = 2 

n = 3 

VMCON 

with Sensitivity 
Derivatives 

1 .033 

1.305 

2.15 

1 .0013 

1 .207 

1 .932 

Max. Potential 
Energy 

with Recurrence 
Procedure 

1 .047 

1.064 

! 

1 .092 

1.005 

i 

1 .024 

1 .048 


Table 2 . Material Distributions for the Optimal Arch Designs of Table 1 Using VMCON 


Type of Loading 

n 

(A e ) 0 PT /(Ae) unif. 

e - 1 

e = 2 

e = 3 

e = 4 

e = 5 

Concentrated 
Load at 
the Crown 

■ 

i 

0.8774 

0.8662 

0.8839 

1 .0452 

1 .340 

2 

0.7036 

0.8370 

0.9963 

1.1577 

1 .3153 

3 

0.5780 

0.8860 

1 .0700 

1.1953 

1.2907 

Uniformly 
Distributed 
Vertical Loading 

■1 

0.9471 

1 .0240 

1 .0455 

1.0080 

0.9760 

2 

0.7662 

0.9526 

1.0663 

1.1080 

1 .1168 

3 

0.6122 | 

0.9285 

1 .0927 

1.1738 

1.2095 


175 































































